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The properties of an innnite-[7 Anderson impurity coupled to both normal and superconducting 
metals is studied using a generalization of the noncrossing approximation which incorporates multiple 
Andreev reflection. Both the cases of a quantum dot and a quantum point contact containing an 
impurity are considered. We find that the magnitude of the Kondo resonance is altered, and that 
structure develops at energies corresponding to the superconducting gap. This leads to observable 
changes in the zero-bias conductance. We also find that magnetic and non-magnetic Kondo effects 
respond in opposite manners to Andreev reflection. 
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I. INTRODUCTION 

Though the Kondo effect is far from being a new prob- 
lem in condensed matter physics, the question of what 
occurs when one places a Kondo impurity on the nor- 
mal side of a normal-superconducting (NS) interface has 
only recently received attention The NS case is 

significantly different from the more studied problem of 
a magnetic impurity in a bulk superconductor, in which 
a lack of low-lying excitations leads to a suppression of 
the Kondo effect. In the present problem, we have both 
the existence of low lying excitations (provided by the 
normal metal) and anomalous pair correlations (a result 
of the proximity effect). The question of what influence 
this combination will have on the dynamically-generated 
Kondo effect is the subject of this work. In particular, 
how will the sharp Kondo resonance at the Fermi energy 
be modified? 

It should be noted that the immediate motivation for 
studying this problem has come from experiment. Recent 
results obtained from semiconductor quantum dots have 
shown clear signs of the Kondo effect ||||, leading to 
the question of what one would expect if a quantum dot 
were now coupled to both normal and superconducting 
leads. In another class of systems, zero-bias conductance 
anomalies seen in metallic quantum point contacts have 
been attributed to scattering off effective non-magnetic, 
two channel Kondo impurities Such point contacts 

can be constructed between normal and superconduct- 
ing metals J9j , providing further relevance for the current 
problem. 



II. THE MODELS 

In this work, we study three different systems: a 
quantum dot coupled to both normal and superconduct- 
ing leads (NS-QDOT), a quantum point contact be- 



tween normal and superconducting metals which con- 
tains a magnetic impurity (M-QPC), and a normal- 
superconducting point contact which contains a non- 
magnetic, two-channel Kondo impurity (NM-QPC). We 
outline in this section the models used for each of these 
devices. 

For the first system, the NS-QDOT, we treat the quan- 
tum dot as an infinite- U Anderson impurity coupled to 
both normal and superconducting leads. The infinite-?/ 
Anderson model has been used successfully in the ab- 
sence of superconductivity to describe the Kondo effect 
in quantum dots; our model represents a natural gener- 
alization to the NS case. Using a slave-boson represen- 
tation, we have: 

H dot ^Ho + e^ flL + W ^ (cl M tff a + /i.e.), (1) 

H o=Y, £ " c lk^ aM + E( Ac L T 4,- fe i + h - c -)- (2) 

a,k,a k 

The cjj t operators here create band electrons, with a 
denoting spin and a = N,S labelling the two leads. A 
represents the pair potential in the superconducting lead. 
The Anderson impurity has bare energy Sd, and is rep- 
resented in the usual manner using auxiliary fermion (/) 
and boson (b) operators; the U — oo constraint of single 
occupancy takes the form + b% = 1. 

For the NS point contact with magnetic impurity sys- 
tem (M-QPC), we use a model consisting of a rectan- 
gular wire of transverse dimensions w, extended in the 
z direction and having a step-function pair potential 
A(z) = Q(z)A. The magnetic impurity is modeled as 
an infinite-?/ Anderson impurity sitting at a point a on 
the NS interface: 

H m = H Q0 + s d J2ftf„ +^^(*l(a)6 t /„ + h.c.) (3) 

a a 

where ^(x) creates a band electron at position x and 
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H Q0 = J dxfe*Ux){-^-E P )* a W + (4) 
A(z)*|(x)*{(x) +hc\ 

For the final system of an NS point contact containing 
a non-magnetic, 2-channel impurity (NM-QPC), we use a 
model similar to that for the M-QPC system. Now, how- 
ever, the impurity represents a two-level system (TLS) 
off which electrons may scatter causing level-flips of the 
TLS fllQfl . In this model, the impurity does not carry or- 
dinary spin, but rather has a pseudo-spin index r which 
can take on one of two values, corresponding to the two 
states of the TLS. The interaction here between impurity 
and conduction electrons is also non-magnetic, meaning 
that electron spin must be conserved. This is accom- 
plished by having the slave-bosons carry a spin index. 
We have: 

H nm = H Q0 + e d J2flL + Wj2^tA a ) b tfr + h.c.) (5) 

where the conduction electron Hamiltonian Hqq is as in 
(H), with the modification that the conduction electron 
operators also carry the pseudo spin index r. We take 
the pairing in (@) to be diagonal in this index, which 
is compatible with the usual interpretation of the the 
pseudo-spin as a parity index. 

It is important to note that there are significant differ- 
ences between the three systems we study, as can be seen 
from their respective Hamiltonians. In the NS-DOT sys- 
tem, the N and S leads are only in contact through the 
impurity; all transport through the system will involve 
it. In the point contact systems, we have the opposite 
situation- the N and the S metals are in perfect contact, 
with the impurity acting only as an additional source of 
scattering at the interface. 

The difference between the two point contact systems 
is also worth emphasizing. In the M-QPC system, our 
impurity is magnetic and thus conduction electron spin 
is not conserved- real spin flips can occur. In the NM- 
QPC system, conduction electron spin is conserved; only 
the auxiliary pseudo-spin index r can undergo spin flips. 
This difference will prove to be crucial, due to the sen- 
sitivity to spin ordering of the singlet pairing in the su- 
perconductor. 

III. EXTENSION OF THE NON-CROSSING 
APPROXIMATION 

For all three of our systems, we calculate the impu- 
rity spectral function (also called the impurity density 
of states) using an extension of the self-consistent non- 
crossing approximation (NCA) p] . The NCA amounts 
to an infinite re-summation of perturbation theory, and 



has been shown to be quantitatively reliable down to tem- 
peratures well below Tk Jll|. We modify the NCA to 
now include multiple- Andreev reflection processes. The 
resulting NS-NCA /-fermion and slave boson propaga- 
tors are given in Fig. 1. The new graphs in our approxi- 
mation are those which contain anomalous propagators; 
they do not appear in the usual NCA. 

Like the original NCA, the NS-NCA is a conserving 
approximation as it may be derived by differentiating a 
generating functional (i.e.. it is <I>-derivable) |ll|. More- 
over, if N is the degeneracy of the Anderson impurity, 
the NS-NCA includes all graphs to order 1/N, including 
new Andreev reflection graphs which first appear at this 
order. As the success of the normal NCA is attributed to 
the fact that it too is exact to order 1/N (in the absence 
of superconductivity), the NS-NCA employed here can 
be viewed as a natural extension to NS Kondo systems 

& 
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FIG. 1. Diagrammatic representation of the NS-NCA. 
Dashed lines are f-fermions, wavy lines are slave bosons, solid 
lines are lead electrons. Double lines indicate a fully dressed 
propagator. Anomalous propagators indicate Andreev reflec- 
tion. 



It should be noted that previous studies of NS Kondo 
systems have either completely neglected changes in im- 
purity dynamics due to superconductivity |ffl], or have 
relied on the equations of motion approach which is for- 
mally only valid at temperatures above the Kondo tem- 
perature In analogy to the usual NCA, we expect 
our approach based on the NS-NCA to be quantitatively 
valid for temperatures well below Tk- 

The Dyson equations pictured diagrammatically in 
Fig. 1 lead to a set of coupled integral equations for 
the /-fermion and slave boson propagators. Letting 
F(uj) = (w-s d - XM)" 1 and B{u) = {w - n(w))- 1 
represent the /-fermion and slave boson retarded propa- 
gators respectively, the equations read: 

E ( W ) = — Ids (p{e)B{u J - s)f(-s) T (6) 
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de'a(e)a(e')B(uj + e)F(w + e + e')B{uj + e') 



n(«) 



2r 



de[p(e)F(Lu + e)f(e)± (7) 



de'a{e)a{e')F(u) + e)B(w + e + e')F(uj + e') 



In these equations, M is the channel degeneracy of the 
impurity- it equals one for both the NS-QDOT and the 
M-QPC systems, and equals two for the NM-QPC system 
(where the two conserved values of electronic spin play 
the role of the two channels). The upper (lower) sign 
corresponds to the M = 1 (AT = 2) systems. p(e) is the 
electronic density of states, T = irW 2 p(0) the bare tun- 
neling rate, / the Fermi distribution function, and a(e) 
is an effective electron-hole coherence parameter defined 
by: 

a(w) =5>n( a K( a )*(M -e„)/(«) (8) 

n 

where u„ and v„ are the usual BCS coherence factors. 




FIG. 2. Diagrammatic argument for the different sign of 
the Andreev reflection graphs in equations (^|) . The first line 
shows a typical double Andreev reflection event for a mag- 
netic impurity. The impurity carries a spin index in this case 
(i.e. the dashed propagator), and we see that both Andreev 
reflections are related by time reversal symmetry. The second 
line shows the situation for a non-magnetic impurity. Here, 
the impurity itself does not carry a spin, but the slave-boson 
(wavy line) does in order to conserve spin. The two Andreev 
reflections in this case are not time reversed partners; the 
second is off by a spin rotation. 



Andreev reflections are not time-reversed partners, but 
are in fact off by a single spin rotation, leading to the 
additional factor of (—1) (recall that the BCS anomalous 
propagator changes sign under a spin rotation). This 
argument is displayed graphically in Fig. |^. A similar 
argument can be made to explain the sign difference oc- 
curring in the slave-boson graph. 

We thus see that the new Andreev reflection processes 
which contribute to the dynamics of the impurity are sen- 
sitive to whether or not the impurity is magnetic; this is 
a direct consequence of the spin-sensitivity of the sin- 
glet pairing in the superconductor. This difference will 
have a significant consequence on how the Kondo effect 
is modified, as will be discussed in what follows. 



IV. RESULTS 



In all three of the systems considered, we choose model 
parameters corresponding to having the Anderson im- 
purity in the Kondo regime, and use a Gaussian with 
half- width D for the normal-state density of states. Our 
choices of e d = -0.67D, T = 0.15D yield T K = 0.0005D 
for a one-channel impurity (i.e. the NS-QDOT and M- 
QPC systems), and Tk = 0.0001D for a two-channel 
impurity (i.e. the NM-QPC system). We have numer- 
ically solved the NS-NCA equations in equilibrium for 
various temperatures and values of the superconducting 
gap. Within the NCA, the impurity spectral function 
Ao{w) can be directly related to the /-fermion and slave 
boson spectral functions: 



AtoM = / ds[e-P £ + e 



}A fa (e)A b (e-cj) (9) 



where the auxiliary particle spectral functions are defined 



by A f = -|lmF, A b 



ilm_B. Note that the equal- 
ity in equation (|9j) reflects a neglect of vertex corrections 
which is consistent with the large- N nature of the NCA. 



It is worth commenting at this point on why the new 
Andreev reflection terms in equations (jfy and (0) enter 
with a different sign for the NM-QPC system compared 
to the NS-QDOT and M-QPC systems. The crucial dif- 
ference has nothing to do with the channel symmetry of 
the impurity (i.e. 2 versus. 1), but rather lies in the 
fact that the impurity-band interaction in the NM-QPC 
is non-magnetic, whereas it is magnetic in the NS-QDOT 
and M-QPC systems. The fact that conduction electron 
spin must be conserved in the non-magnetic case has a 
direct consequence on the sign of the Andreev-reflection 
graph appearing in the /-fermion self energy- the two 



A. NS Quantum Dot 



In Fig. ||, we plot calculated spectral functions for the 
NS-QDOT system for several values of the superconduct- 
ing gap A. In the case of no superconductivity (A = 0), 
we see as expected the emergence of the sharp Kondo 
resonance at the Fermi energy. As the superconductivity 
is gradually turned on, several interesting modifications 
of this resonance are observed. 



3 




-20 -15 -10 -5 5 10 15 20 

Energy(T K ) 

FIG. 3. NS QDOT spectral function A da (u) for various 
values of A at T = 0.5T K 



The Kondo resonance does not vanish completely as A 
increases, though the spectral weight associated with the 
resonance clearly decreases, indicating a partial suppres- 
sion of the Kondo effect. This is to be expected when one 
recalls that the ground state of a magnetic, one-channel 
Kondo impurity is one in which the local impurity screens 
its spin with low-energy conduction electrons. Turning 
on the superconductivity leads to a gap in the supercon- 
ducting lead, and thus there are fewer low energy ex- 
citations available to the impurity to use in screening. 
It thus becomes more difficult to form the spin-screened 
Kondo ground state, resulting in a diminishing of the 
Kondo effect. A similar result for the NS QDOT is found 
using slave-boson mean- field theory Q, where the effec- 
tive Kondo temperature is found to fall as A is increased 
above Tk- 

Also of interest are the sub-peaks which develop in 
the impurity spectral function at roughly ±A. These 
peaks are a result of the Kondo effect, and do not appear 
at higher temperatures. Recalling that the the impu- 
rity spectral function at temperature can be viewed as 
a local particle addition/removal spectrum, these peaks 
indicate that one can create a superconducting-like ex- 
cited state with excitation energy A by adding a particle 
to the quantum dot. This indicates that correlations are 
developing between the superconductor and the impurity, 
implying that the superconducting electrons do indeed 
participate in the Kondo effect, even for A >> Tk- This 
can be substantiated by using the large- -/V variational ap- 
proach to calculate the approximate ground state of 
the NS-QDOT. One finds that there are superconduct- 
ing quasiparticles present in the interacting ground state, 
with a weight that scales as Tr-/A for large A. 

As we find that a Kondo resonance persists in the NS- 
QDOT system, a natural question to ask is whether there 
is any resulting enhancement of the zero-bias conduc- 
tance of the device, as is seen in the absence of supercon- 



ductivity. The equilibrium spectral functions obtained 
for the NS-QDOT can be used to calculate the zero-bias 
Andreev conductance of the device, using an approximate 
formula derived in 0] : 

G NS = / ^Re£f 2 Im(Df 2 D£)[-|-] (10) 

Here, D^ A ^ labels a component of the retarded (ad- 
vanced) impurity green function in Nambu-space, and 
£^" 4 ' ) is the corresponding self-energy. We have com- 
puted © within the NS-NCA, evaluating the anoma- 
lous impurity green function from the diagram in Fig. 3 
without use of the standard "elastic" approximation p9|. 




FIG. 4. Diagram used to calculate anomalous impu- 
rity propagator within the NCA. Double-dashed lines rep- 
resent fully dressed f-fermion, double wavy lines represent 
dressed slave boson propagators; both are calculated within 
the NS-NCA. The solid line represents an anomalous lead 
electron propagator. 

For Tk < A < T and temperatures as low as one- 
half T K , we find that the NS QDOT does not exhibit 
a Kondo-induced enhancement of the zero-bias conduc- 
tance, despite the presence of a Kondo resonance. This 
can be understood by the fact that in the infinite U limit , 
there is a suppression of pairing at the dot compared to 
the U = case. For energies and voltages smaller than 
A, the only processes to contribute to the current of the 
NS QDOT are those in which two particles tunnel coher- 
ently through the dot and enter the superconductor as 
a pair. The suppression of pairing at the dot thus sup- 
presses the zero bias conductance, despite the fact that 
there still is a Kondo resonance in the impurity spectral 
function. Note that the same conclusion is obtained using 
a correct equations-of-motion approach (see the erratum 
in Ref. @). 

B. NS Point Contact Systems 

We proceed to discuss results for the NS point contact 
with magnetic impurity (M-QPC) and NS point contact 
with non-magnetic impurity (NM-QPC) systems. Shown 
in Fig. |^ is a plot of the impurity spectral function for 
the NM-QPC system for various values of A and at a 
temperature below the bare Kondo temperature; Fig. ^ 
shows a similar plot for the M-QPC system. 
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Several features are noteworthy in these plots. First, 
we notice that in both systems the introduction of su- 
perconductivity does not cause a significant change in 
the amount of spectral weight in the Kondo resonance. 
This indicates that there is no major suppression of the 
Kondo effect here, in contrast to the behaviour observed 
in the NS-QDOT. This is to be expected- for a clean, 
ballistic NS point contact, an energy gap does not form 
at the position of the impurity (which is at the interface) 
p5[ , and consequently there is no shortage of the low 
energy excitations needed to generate the Kondo effect. 
This differs from the NS QDOT, where the impurity is 
coupled directly to a bulk superconductor having a fully 
formed gap. 




Energy(T K ) 



FIG. 5. Impurity spectral function Ad a (uj) for a 

non-magnetic two channel Anderson impurity in an NS QPC 
for various values of A at T — 0.66Tk- Note that the intro- 
duction of superconductivity reduces the height of the Kondo 
resonance. 
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FIG. 6. Impurity spectral function Ada (w) for a magnetic 
one channel Anderson impurity in an NS QPC at T = 0.50Ta-. 
Note here that the introduction of superconductivity increases 
the height of the Kondo resonance. 

More interestingly, we find that in both point contact 
systems there is a small modification of the height of the 
Kondo peak when A is increased from 0. In the NM- 
QPC, we find that the Kondo peak is reduced, whereas 
the opposite is observed in the M-QPC- here, the Kondo 
peak is slightly enhanced. 

This behaviour can be understood heuristically as fol- 
lows. First, note that by introducing superconductivity 
into the system, we introduce new processes which fur- 
ther couple the impurity to the conduction electrons (i.e. 
the multiple Andreev reflection processes which corre- 
spond to the new graphs included in the NS-NCA) . These 
additional processes modify the effective coupling T be- 
tween the impurity and conduction electrons; an increase 
of this coupling will enhance the Kondo effect, while a 
reduction will suppress it |l6[ . In the present case, the 
sign of this modification depends on whether or not the 
impurity is magnetic- as discussed in Sec. Ill, the new 
Andreev-reflection graphs enter with opposite signs in 
these two cases due to the sensitivity of the Cooper pair- 
ing to spin ordering. This leads directly to the opposite 
behaviour of the Kondo peak in the two cases. 

The equilibrium impurity spectral functions can be 
used to derive the zero bias conductance for the point 
contact systems. We compute the current using the non- 
equilibrium Keldysh technique as was done in fl7],[ll| , re- 
stricting ourselves to the regime where V, T << A. Using 
particle-hole symmetry, the zero bias conductance takes 
the form: 

4e 2 f df 
G NS = 2G N - —nT I dojA d {io)[-^-{uj)] (11) 

where Gat is the normal-state Sharvin conductance of 
the point contact in the absence of the impurity. Details 
of this calculation are presented in the appendix. Note 
that we have neglected a term involving the anomalous 
impurity spectral function which is for T = and is 
in general much smaller than the Ad{u>) term due to the 
strong suppression of on-site pairing. 

The first term in (^) expresses the usual doubling of 
the conductance expected for an ideal NS point contact, 
while the second term represents the correction due to 
the impurity (SGns)- As in the N case, the formation of 
the Kondo resonance will lead to a small suppression of 
the point contact conductance around zero voltage. The 
correction term due to the impurity in the absence of su- 
perconductivity (SGn) has the same form as the second 
term in (|lT|), but with an additional factor of 1/2 Jl8| . 
If Andreev reflection did not effect the dynamics of the 
impurity, Ad(to) would be identical in both the NS and 
N cases. We would thus expect the magnitude of the NS 
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conductance anomaly to be twice that of the N anomaly. 

In Fig. ^ we plot the parameter g = (SGns — 
2SGn)/SGn as a measure of the size of the NS QPC 
zero bias anomaly. We see that as T is lowered below 
Tk, g becomes non-zero in both the cases of a magnetic 
and a non-magnetic impurity indicating that Andreev 
reflection docs indeed have an observable impact on the 
Kondo effect in these systems. Note also that the sign 
of g in both cases is in accord with our discussion of the 
respective impurity spectral functions- in the magnetic 
case, the NS zero bias anomaly is more than two times 
the size of the anomaly in the normal state, while in 
the non-magnetic case, the NS anomaly is smaller than 
twice the normal state anomaly. Again, magnetic and 
non-magnetic Kondo effects respond in opposite ways to 
Andreev reflection. Significant too is the fact that the 
effect is large enough (±20% at the lowest temperatures 
tested) that it should be experimentally accessible. This 
behaviour could provide another test of the two-channel 
Kondo explanation of zero bias conductance anomalies 
seen in metal point contacts. 



Magnetic 1 -Channel. A = 25T K 
Non-magnetic 2-channel, A = 25T K 
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FIG. 7. Parameter g = (SGns - 2SGn)/5Gn for the NS 
QPC as a function of temperature. The non-zero value of 
g is a consequence of changes in impurity dynamics due to 
Andreev reflection. 



V. CONCLUSION 

We have studied the properties of three NS Kondo sys- 
tems using an extension of the NCA which incorporates 
Andreev reflection. We find the Kondo effect is indeed 
modified by superconductivity in each system, leading to 
changes both in the impurity spectral function and in the 
zero-bias conductance. In the case of the NS-QDOT, su- 
perconductivity causes an overall decrease in the spectral 
weight of the Kondo resonance, and leads to the forma- 
tion of new Kondo sub-peaks at roughly ±A. There is no 
resulting enhancement of the zero-bias Andreev conduc- 
tance due to the suppression of on-site pair correlations 
in the large-?/ limit. In the point contact with impu- 



rity systems, we find that superconductivity enhances the 
Kondo peak height in the case of a magnetic impurity, 
and suppresses it in the case of a non-magnetic impurity. 
We attribute this difference to the sensitivity of Cooper 
pairing to spin ordering. 
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APPENDIX A: CONDUCTANCE OF AN NS 
POINT CONTACT WITH IMPURITY 

We derive in this appendix Eq. (^) for the conduc- 
tance of a constriction between a normal metal and a 
superconductor which contains an interacting impurity. 
Our technique is similar to that used in jl8| to calculate 
the current of a normal metal point contact with an im- 
purity. Here, however, the "trick" of combining currents 
calculated to the left and right of the constriction to elim- 
inate the impurity lesser function cannot be used- we do 
not wish to deal with the complexities of calculating the 
supercurrent on the S side of our constriction. We in- 
stead focus solely on the current in the normal region, 
and use an approximate electron-hole symmetry which is 
present at low energies to eliminate the impurity lesser 
function. 

We begin by writing the current in terms of the lesser 
Green function: 



/ = -~2ieh 



djj 
2^ 



d x i 



d z - d z 
2mi 



G < 



(Al) 



where G < (x, x';w) is the Fourier transform in time of: 



G < (x,t;x',0) =i/*}(x',0)* T (x,i) 



(A2) 



Note that we use only the 11 component of the matrix 
Nambu Green function in computing the current; a fac- 



tor of two is included in (Al) to account for spin. Note 
also that in what follows we neglect the pseudo-spin in- 
dex r which appears in the case of a two-channel Kondo 
impurity; our discussion may nevertheless be applied to 
this case by simply averaging over this index. 

The next step in the derivation is to express the current 
in terms of the quasi-particle operators Jnvtfna which 
appear when one makes a Bogolubov-de Gennes (B-dG) 
transformation to diagonalize the non-impurity parts of 
the Hamiltonian (i.e. *|(x) = w «( x )7r;T - <( x blx 
etc.). This leads to the following identification: 



G<(x,x'; 



with: 



0i,m(x) G<„(w) 



(A3) 
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6i, m (x) = 

^2,m(x) = 



«m(x) 

%(x) 



^ >T (0)7m,T(*)^ (7na(0)7m,T(*)> 

jit (0)7,1,1 W) (7^(0)7^4(0' 



(A4) 
(A5) 

(A6) 



In Eq. (A3), i and j are indices in particle-hole space, 
whereas m and n label different quasiparticle modes; we 
use the* symbol to denote structure in particle-hole space. 
All repeated indices are to be summed over. Note that 
the B-dG wavefunctions u and v have a particularly sim- 
ple form for mode energies e less than the gap A [f20| . 
Letting rh index the transverse wavefunctions ipm{x-±), 
we have for incident-electron modes: 



h,m — 01,(m,e) 



exp(ifc e (e)z) 
a(e) eyLp(ik h {e)z) 



^m(Xj.) (A7) 



where 



(Ep + e — E m ), 



kt, 



2m 1 
IP 



(Ep — e—E m ) and a — exp(— i arccos(e/A)) is 
the Andreev reflection phase. Ep denotes the Fermi en- 
ergy, while E m is the energy of the transverse mode rh. 
Incident-hole wavefunctions have a similar form: 



? l,(m,e) 



a(e) exp(—ik e (e)z) 



fe(xi) (A8) 



The problem now becomes one of computing the mode- 
space lesser functions G^\ n (oj). As the only non-trivial 
interaction in our system is at the impurity itself (the 
quasiparticles only interact through a hopping term), we 
may express the exact self-energy appearing in the Dyson 
equation for G< n (o;) in terms of D(ui), the Nambu ma- 
trix Green function describing the impurity. 

The Dyson equation takes the form: 



C < = 5 
^ mn 



a 



mn 
■0,R 



G 



0,< 



r << 4. 



Q.A 



(A9) 



G 



0,< s^A 



y A g 

mn ' 



Q.A 



where the superscript indicates a Green function of the 
non-interacting system, and the R (A) superscript in- 
dicates a retarded (advanced) Green function. All func- 
tions in the above formula are to be evaluated at the same 



frequency. The self-energy matrix appearing in (AE) is 
simply: 



[<#a,m( a ) 



,n(a) 

J 3 

(A10) 



with similar formulas for the advanced and lesser self- 
energies. Recall that a indicates the position of the im- 
purity. This simple form for the self-energy reflects the 



fact that quasiparticles only interact by hopping on and 
off the impurity site. 



We proceed by substituting Eq. (A9) into Eq 



We evaluate the resulting expression in the limit where 
z, z' both tend to —00 (i.e. deep in the normal region). 
As the wave functions u m (x) and V m (x) have a plane 
wave form, the terms required to calculate the current are 
greatly simplified. A further simplification arises when 
we restrict attention only to those quasiparticle modes 
with energies less than A; this is valid for voltages and 
temperatures which are much smaller than A. 

After a fair amount of algebra, we obtain a relatively 
simple formula for the current. Ignoring terms involving 
the anomalous impurity Green functions for the moment 
and writing I = Iq + 51, we have: 



fdu; £ (^H - //.M) ( A11 ) 



77 ' lL ' 



51 



2e 



+ ___£>< ( w ) - A 22 (u)f h (uj)) (A12) 
2m J 



In the above equations, f e {^) is the distribution func- 
tion describing electrons, while fh(u>) = 1 — f e (—u>) de- 
scribes holes. Aij(oj) = ^Im(Df!j (a;)) is an impurity spec- 
tral function. T(u) = ttW 2 N(lo) is an energy-dependent 
tunneling rate, with N(ui) being the density of states. 
Of course, the energy integrals in both these equations 
should be restricted to values less than A due to the ap- 
proximations that have been made; this will occur natu- 
rally for sufficiently small temperatures and voltages. 



Io in Eq. ( All ) is the point contact current in the 
absence of the impurity. Taking f e {u) to be the shifted 
Fermi function f(u — eV), we find that Go = ^7 = 

^r-N m , where N m is the number of transverse modes 
with energies below the Fermi energy. This is precisely 
twice the usual Sharvin conductance, and represents the 
current-doubling effect of Andreev reflection j2Q] . 

Eq. (|A12| ) for the modification of the current due to 
the impurity is still in a somewhat unwieldy form, due to 
the appearance of the impurity lesser functions. While 
the NS-NCA technique discussed in this paper allows the 
calculation of spectral functions, it does not allow one to 
calculate the distribution functions which are necessary 
to compute the lesser Green function. A similar problem 
occurs when considering an impurity in a normal point 
contact fll8| , or more generally, the current through an 
interacting region p7| . In these cases, one can eliminate 
the lesser function by taking a suitable linear combina- 
tion of the current computed on the left of the interact- 
ing region and that computed on the right. Here, no 
such simplification is possible, as computing the current 
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within the superconductor is difficult due to the presence 
of supercurrents. 

A means for eliminating the lesser function does how- 
ever exist. We note that the electron impurity Green 
functions Dn(w) are related to the hole impurity Green 
functions Z^M- It follows easily from the definitions of 
the functions that: 



D^(uj) = -D> (-w) 
A 22 (u) = An(-uj) 



(A13) 
(A14) 



where the greater Green function is defined by: 



D> 1 (t)=-i(d r (t)d\(0) 



(A15) 



We use Eqs. ( |A13| ) and ( |A14|) to substitute for the 
third and forth terms in Eq. ( |A12 ) for SI. Further, we 
change the integration variable for these terms from w 
to — oj, and make the assumption that T(u>) is an even 
function. The latter is not too severe as long as the tem- 
perature and voltage are sufficiently small, as this ensures 
a small range of u> integration. Finally, we use the rela- 
tions fh(—(j) = 1 — / e M and Dj^ — = —2mAu to 
obtain: 



SI 



4e 



du)(7rT(u/))Ad(u>) 



fe{oj) - fh(u) 



(A16) 



Note that Ad in the above formula is the same as An. 
Taking T to be a constant and differentiating with re- 
spect to voltage at V=0 yields the formula ( pd| ) for the 
impurity contribution to the point contact zero-bias con- 
ductance. 

We now comment on terms in the current which involve 
anomalous impurity green functions. Retaining these 
terms and performing manipulations similar to those de- 
scribed above leads to a second correction term to the 
point contact current: 

xt 4e [a 1 rv a < \\ f/eM _ AM 

Slanorn = J~ / dui(nT(u}))( — Ai2(u))) I 

(A17) 

The change to the zero-bias conductance arising from 
this term is at temperature, and is generally much 



smaller than SI in (A16) for temperatures and voltages 



smaller than A. In the case of a strong-!/ Anderson im- 
purity, Ai2(lo) is strongly suppressed due to the suppres- 
sion of on-site pair correlations. Even in the absence of an 
on-site repulsion, ^12(0;) is small for u> < A, as there are 
no non-evanescent quasiparticles in the gap. The pres- 
ence of the factor further reduces the contribution of 
this term for small voltages and temperatures. We thus 
neglect the contribution of this term in computing the 



point contact zero-bias conductance. Note that we have 
evaluated numerically the significance of this term in a 
few test cases by calculating A12 within the NS-NCA; it 
typically changes the magnitude of the zero-bias anomaly 
by less than 0.5%. 
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